1 



study of J/il^ ujpp at BESIII 



2 M. Ablikim\ M. N. Achasov*^, O. Albayrak^ D. J. Ambrose^^, F. F. An\ Q. An-*", J. Z. Bai\ R. Baldini Ferroli"-*, 

3 Y. Ban2^ J. Becker^ J. V. Bennettl^ M. Bertani"-*, J. M. Bian^^ E. Boger^'''", O. Bondarenko^^ , I. Boyko^^ R. A. Briere^ 

4 V. Bytev", H. Cai**, X. Cai\ O. Cakir^'^'*, A. Calcaterra^^-^, G. F. Cao\ S. A. Cetin=*^^, J. F. Chang\ G. Ghelkov"'", 

5 G. Chen\ H. S. Chen\ J. C. Chen\ M. L. Chen\ S. J. Chen^", X. Chen^'', Y. Chen^ Y. B. Clien\ H. P. Cheng", 

6 Y. R Chu\ D. Cronin-Hennessy='*, H. L. Dai\ J. P. Dai^ D. Dedovich^^ Z. Y. Deng\ A. Denig^^ I. Denysenko"■^ 

7 M. Destefanis*^-*'"^'^, W. M. Ding^*, Y. Ding^^, L. Y. Dong\ M. Y. DongS S. X. Du''^ J. Fang\ S. S. Fang\ L. Fava*^^'"^'' 

8 C. Q. Feng''^ P. Friedel^, C. D. Fu\ J. L. Fu^"*, O. Fuks^^'", Y. Gao^^ C. Geng''°, K. Goetzen^ W. X. Gong\ W. Gradl^**, 

9 M. Greco*^'*'-*^^, M. H. Gu\ Y. T. Gu^ Y. H. Guan^^ A. Q. Guo=^ L. B. Guo^^, T. Guo^^ Y. P. Guo^^, Y. L. Han\ 

10 F. A. Harris^^ K. L. He\ M. Re\ Z. Y. He2^ T. Held^ Y. K. Heng\ Z. L. Hou\ C. Hu^^, H. M. Hu\ J. F. Hu■'^ T. Hu^ 

11 G. M. Huang*, G. S. Huang*", J. S. Huang^^ L. Huang\ X. T. Huang2^ Y. Huang^*, Y. P. Huang\ T. Hussain-*^ C. S. Ji*° 
- - - 12 Q. Ji\ Q. P. JF', X. B. Ji\ X. L. Ji\ L. L. Jiang\ X. S. Jiang\ J. B. Jiao2^ Z. Jiao", D. P. Jin\ S. Jin\ F. F. Jing='^ 
^ i 13 N. Kalantar-Nayestanaki^", M. Kavatsyuk^", B. Kopf^ M. Kornicer^'', W. Kuehn^^ W. Lai\ J. S. Lange='^ P. Larin", 
^ . 14 M. Leyhe^ C. H. U\ Cheng Li*°, Cui Li*", D. M. Li''^ F. Li\ G. Li\ H. B. U\ J. C. Li\ K. Li^°, Lei Li\ Q. J. Li\ 
^ ' 15 S. L. Li\ W. D. Li^ W. G. Li\ X. L. Li^^, X. N. Li\ X. Q. LP, X. R. Li^^ Z. B. LP, H. Liang"", Y. F. Liang^°, 
^ ; 16 Y. T. Liang■^^ G. R. Liao^^ X. T. Liao\ D. Lin", B. J. Liu\ C. L. Liu^ C. X. LiuS F. H. Liu^^ Fang LiuS Feng Liu*, 

>>. 17 H. Liu\ H. B. Liu^, H. H. Liu", H. M. Liu\ H. W. Liu\ J. P. Liu**, K. Liu^^ K. Y. Liu^^ Kai Liu^^ P. L. Liu^^ Q. Liu^*^ 
C3 ■ 18 S. B. Liu*", X. Liu2\ Y. B. Livi^^ , Z. A. Liu*, Zhiqiang Liu\ Zhiqing Liu*, H. Loehner^", G. R. Lu*^ H. J. Lu**, J. G. Lu*, 
^ ; 19 Q. W. Lu2^ X. R. Lu^^ Y. P. Lu\ C. L. Lno^^, M. X. Luo*^ T. Luo^^ X. L. Luo\ M. Lv*, C. L. Ma='^ F. C. Ua^^ , 
. 20 H. L. Ma\ Q. M. Ma*, S. Ma\ T. Ma\ X. Y. Ma\ F. E. Maas**, M. Maggiora*^'*'*^'^ , Q. A. MaIik*^ Y. J. Mao^^ 
CN ■ 21 Z. P. Mao\ J. G. Messchendorp^", J. Min*, T. J. Min*, R. E. Mitchell*^ X. H. Mo*, H. Moeini^", C. Morales Morales**, 
CN ' 22 K. Moriya*^ N. Yu. Muchnoi^ H. Muramatsu^^ Y. Nefedov*^ C. Nicholson^^ I. B. Nikolaev'^, Z. Ning*, S. L. Olsen^^ 
! 23 Q. Ouyang*, S. Pacetti*'^'^ , J. W. Park^^ M. Pelizaeus^ H. P. Peng*°, K. Peters^ J. L. Ping^^, R. G. PingS R. Poling3^ 
'77" 24 E. Prencipe", M. Qi^*, S. Qian\ G. F. Qiao^^ L. Q. Qin^^, X. S. Qin*, Y. Qin^^, Z. H. Qin*, J. F. Qiu*, K. H. Rashid*^ 
5^ ■ 25 G. Rong*, X. D. Ruan^ A. Sarantsev*^'^ B. D. Schaefer*®, M. Shao*°, G. P. Shen^'^-'', X. Y. Shen*, H. Y. Sheng*, 
I ; 26 M. R. Shepherd*^ W. M. Song*, X. Y. Song*, S. Spataro*^'*'*^'^ , B. Spruck^^ D. H. Sun\ G. X. Sun*, J. F. Sun^^ S. S. Sun* 
Ph- 27 Y. J. Sun*°, Y. Z. Sun*, Z. J. Sun*, Z. T. Sun*", C. J. Tang^", X. Tang*, I. Tapan^*^, E. H. Thorndike='^ D. Toth=**, 
<D ■ 28 M. Ullrich^^ I. Uman^*^, G. S. Varner'^^ B. Q. Wang2^ D. Wang2^ D. Y. Wang2^ K. Wang*, L. L. WangS L. S. Wang*, 
|-C ; 29 M. Wang2^ P. Wang*, P. L. WangS Q. J. Wang*, S. G. Wang2^ X. F. Wang3^ X. L. Wang*", Y. D. Wang*''-'', Y. F. Wang* 
. 30 Y. Q. Wang■*^ Z. Wang*, Z. G. Wang*, Z. Y. Wang*, D. H. Wei^ J. B. Wei^^ , P. WeidenkafT**, Q. G. Wen*", S. P. Wen*, 
' 31 M. Werner^^ U. Wiedner^ L. H. Wu*, N. Wu*, S. X. Wu*", W. Wu^^, Z. Wu*, L. G. Xia^^ Y. X Xia*^ Z. J. Xiao2^ 
^ ; 32 Y. G. Xie*, Q. L. Xiu*, G. F. Xu*, G. M. Xu^^ Q. J. Xu*", Q. N. Xu^^ X. P. Xu^*, Z. R. Xu*", F. Xue*, Z. Xue*, L. Yan*" 
QO ! 33 W. B. Yan*", Y. H. Yan*^ H. X. Yang*, Y. Yang*, Y. X. Yang^ H. Ye*, M. Ye*, M. H. Ye^ B. X. Yu*, C. X. Yu2^ 
Q ■ 34 H. W. Yu2^ J. S. Yu^*, S. P. Yu2^ C. Z. Yuan*, Y. Yuan*, A. A. Zafar*^ A. Zallo*''^, S. L. Zang^*, Y. Zeng*^ B. X. Zhang* 
^ ■ 35 B. Y. Zhang*, G. Zhang^*, G. G. Zhang\ D. H. Zhang*, H. H. Zhang^^ H. Y. Zhang*, J. Q. Zhang*, J. W. Zhang*, 
m ; 36 J. Y. Zhang*, J. Z. Zhang*, LiLi Zhang*^ R. Zhang^^ S. H. Zhang*, X. J. Zhang*, X. Y. Zhang^*, Y. Zhang*, Y. H. Zhang* 
• ■ 37 Z. P. Zhang*", Z. Y. Zhang**, Zhenghao Zhang*, G. Zhao*, H. S. Zhao*, J. W. Zhao*, K. X. Zhao'^^ , Lei Zhao*", Ling Zhao* 
W ■ 38 M. G. Zhao2^ Q. Zhao*, S. J. Zhao*®, T. C. Zhao\ X. H. Zhao^*, Y. B. Zhao*, Z. G. Zhao*", A. Zhemchugov^^'", B. Zheng** 
39 J. P. ZhengS Y. H. Zheng^®, B. Zhong2^ L. Zhou*, X. Zhou**, X. K. Zhou3^ X. R. Zhou*", G. Zhu*, K. Zhu*, K. J. Zhu*, 
LJ ' 40 S. H. Zhu*, X. L. Zhu3^ Y. G. Zhu*", Y. M. Zhu^^ Y. S. Zhu*, Z. A. Zhu*, J. Zhuang*, B. S. Zou*, J. H. Zou* 



X 



(BESIII CoUaboration) 
^ Institute of High Energy Physics, Beijing 100049, People's Republic of China 



2 



Bochum Ruhr- University, 0-44780 Bochum, Germany 

3 



, 44 Carnegie Mellon University, Pittsburgh, Pennsylvania 15213, USA 

" " 45 * Central China Normal University, Wuhan 430079, People's Republic of China 

46 ^ China Center of Advanced Science and Technology, Beijing 100190, People's Republic of China 

47 ® G.I. Budker Institute of Nuclear Physics SB RAS (BINP), Novosibirsk 630090, Russia 

48 ' GSI Helmholtzcentre for Heavy Ion Research GmbH, D-64291 Darmstadt, Germany 

49 * Guangxi Normal University, Guilin 54IOO4, People's Republic of China 

50 " GuangXi University, Manning 530004, People's Republic of China 

51 ^" Hangzhou Normal University, Hangzhou 310036, People's Republic of China 

52 ** Helmholtz Institute Mainz, J ohann- Joachim- Becher-W eg 45, D-55099 Mamz, Germany 

53 *^ Henan Normal University, Xinxiang 453007, People's Republic of China 

54 *^ Henan University of Science and Technology, Luoyang 471003, People's Republic of China 

55 ** Huangshan College, Huangshan 245000, People's Republic of China 

56 *^ Hunan University, Changsha 410082, People's Republic of China 

57 *® Indiana University, Bloommgton, Indiana 47405, USA 

58 *'' (A)INFN Laboratori Nazionali di Frascati, 1-00044, Frascati, 

59 Italy; (B)INFN and University of Perugia, 1-06100, Perugia, Italy 

60 ** Johannes Gutenberg University of Mainz, J ohann- Joachim- Becher-W eg 45, D-55099 Mainz, Germany 

61 *" Joint Institute for Nuclear Research, 141980 Dubna, Moscow region, Russia 



2 



62 KVI, University of Groningen, NL-9747 A A Groningen, The Netherlands 

63 Lanzhou University, Lanzhou 730000, People's Republic of China 

64 Liaoning University, Shenyang 110036, People's Republic of China 

65 Nanjing Normal University, Nanjing 210023, People's Republic of China 

66 Nanjing University, Nanjing 210093, People's Republic of China 

67 Nankai University, Tianjin 300071, People's Republic of China 

68 ■^^ Peking University, Beijing 100871, People's Republic of China 

69 Seoul National University, Seoul, 151-747 Korea 

70 Shandong University, Jinan 250100, People's Republic of China 

71 Shanxi University, Taiyuan 030006, People's Republic of China 

72 Sichuan University, Chengdu 610064, People's Republic of China 

73 Soochow University, Suzhou 215006, People's Republic of China 

74 Sun Yat-Sen University, Guangzhou 510275, People's Republic of China 

75 Tsinghua University, Beijing 100084, People's Republic of China 

76 (A) Ankara University, Dogol Caddesi, 06100 Tandogan, Ankara, Turkey; (B)Dogus 

77 University, 34722 Istanbul, Turkey; (C)Uludag University, 16059 Bursa, Turkey 

78 Universitaet Ciessen, D-35392 Giessen, Germany 

79 University of Chinese Academy of Sciences, Beijing 100049, People's Republic of China 

80 University of Hawaii, Honolulu, Hawaii 96822, USA 

81 ^® University of Minnesota, Minneapolis, Minnesota 55455, USA 

82 '^^ University of Rochester, Rochester, New York 14627, USA 

83 University of Science and Technology of China, Hefei 230026, People's Republic of China 

84 University of South China, Hengyang 421001, People's Republic of China 

85 University of the Punjab, Lahore- 54590, Pakistan 

86 (A) University of Turin, 1-10125, Turin, Raly; (B)University of Eastern 

87 Piedmont, 1-15121, Alessandria, Raly; (C)INFN, 1-10125, Turin, Italy 

88 Wuhan University, Wuhan 430072, People's Republic of China 

89 Zhejiang University, Hangzhou 310027, People's Republic of China 

90 Zhengzhou University, Zhengzhou 450001, People's Republic of China 

91 " Also at the Moscow Institute of Physics and Technology, Moscow 141700, Russia 

92 On leave from the Bogolyubov Institute for Theoretical Physics, Kiev 03680, Ukraine 

93 " Also at the PNPI, Gatchina 188300, Russia 

94 Present address: Nagoya University, Nagoya 464-8601, Japan 

95 (Dated: May 23, 2013) 



The decay J/ip uipp has been studied, using 225.3 x 10® J/ip events accumulated at BESIII. 
No significant enhancement near the pp invariant-mass threshold (denoted as X{pp)) is observed. 
The upper limit of the branching fraction B{J /ip — ?■ u}X(pp) ^pp) is determined to be 3.9 x 10~® 
at the 95% confidence level. The branching fraction of J/tp — > ujpp is measured to be B{J/4> — > 
u)pp) = (9.0 ± 0.2 (stat.) ± 0.9 (syst.)) x 10"''. 

96 PACS numbers: 13.25.Gv, 12.39.Mk, 13.75.Cs 
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An anomalous enhancement near the threshold of the 
pp system, namely X{pp), was first observed by the BE- 
SII experiment in the radiative decay J/ip ^pp 
and it was recently confirmed by the CLEO and BESIII 
experiments 0-111. In the BESII experiment, its mass IS 



measured to be 1859Iio (stat.)l25 (syst.) MeV/c^ and 
the total width is T < SOMcV/c^ at the 90% confi- 
dence level (C.L.). While in the BESIII experiment, a 
partial wave analysis (PWA) with a correction for the 
final-state interaction (FSI) is performed, and the spin- 
parity of X{pp) is determined to be its mass is 
183215^ (stat.)ti7 (syst.) MeV/c^ and the total width 
is r < 76MeV/c2 at the 90% C.L. [|. 

The discovery of X{pp) stimulated a number of the- 
oretical interpretations and experimental studies [5l-[l6|. 
There is no experimental evidence of such an enhance- 



ment in other quarkonium decays, e.g. J/ip — s- ir^pp [l| 
or T{2S) jpp Q. In t/j{2S) jpp, the recent BESIII 
measurement shows a relative production rate to that of 
J/ij decays of i? = 5.08% Q. A number of theoretical 
speculations have been proposed to interpret the nature 
of this structure, including baryonium [9|- 1 ll . a multi- 
quark state [m or mainly a pure FSI |13. Il4l|. It was 
proposed to associate this enhancement with a broad en- 
hancement observed in B meson decays [U [3 or a new 
resonance X(1835) in J/-0 — ^ jTT~^n~r] decay at BE- 
SII [11. 

The investigation of the near-threshold pp invariant 
mass spectrum in other J/ip decay modes will be helpful 
in understanding the nature of the observed structure. 
The decay J/ip — >■ wpp restricts the isospin of the pp 
system, and it is helpful to clarify the role of the pp 
FSI. The BESII collaboration studied J/^p ^ Lupp via lu 
decaying to tt'^tt+tt" with a data sample of 5.8 x 10^ J/tp 
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132 events [6[ . No significant signal near the threshold of the 

133 pp invariant-mass spectrum was observed and an upper 

134 limit on the branching fraction of J/i/j — > ijjX{pp) — > Lopp 

135 was determined to be 1.5 x 10^^ at the 90% C.L., which 
135 disfavored the interpretation of a pure FSI effect giving 

137 rise to the X{pp). In this paper, the analysis oi J/^p — > 

138 upp via the decay channel uj — > 77r'^ is presented, based 

139 on a data sample of (225.3 ± 2.8) x 10^ J/i) events [13 

140 accumulated with the BESIII detector. Searching for the 

141 X{pp) in the decay mode J/il) ^ Lopp jir'^pp has a 

142 particular advantage: a low irreducible background from 

143 A'^* is expected. The channel J/^p — > cupp — > mnr'^pp has 

144 irreducible background from various A'^* decays and A 

145 decays, where interferences may have a large impact on 

146 the uncertainty of the measurements. 

147 BESIII/BEPCII (U is a major upgrade of the BESII 

148 experiment at the BEPC accelerator for studies of 

149 hadron spectroscopy and r-charm physics The de- 

150 sign peak luminosity of the double-ring e^e~ collider, 

151 BEPCII, is 10^3 cm-^s-i at beam currents of 0.93 A. The 

152 BESIII detector with a geometrical acceptance of 93% 

153 of 47r, consists of the following main components: 1) a 

154 small-celled, helium-based main drift chamber (MDC) 

155 with 43 layers. The average single wire resolution is 

156 135 //m, and the momentum resolution for IGeV/c^ 

157 charged particles in a IT magnetic field is 0.5%; 2) 

158 an electromagnetic calorimeter (EMC) made of 6240 

159 Csl (Tl) crystals arranged in a cylindrical shape (barrel) 

160 plus two end-caps. For 1.0 GeV photons, the energy res- 

161 olution is 2.5% in the barrel and 5% in the end-caps, and 

162 the position resolution is 6 mm in the barrel and 9 mm in 

163 the end-caps; 3) a Time-Of-Flight system (TOE) for par- 

164 tide identification (PID) composed of a barrel part made 

165 of two layers with 88 pieces of 5 cm thick, 2.4 m long plas- 

166 tic scintillators in each layer, and two end-caps with 48 

167 fan-shaped, 5 cm thick, plastic scintillators in each end- 

168 cap. The time resolution is 80 ps in the barrel, and 110 

169 ps in the end-caps, corresponding to a K/tt separation 

170 by more than 2a for momenta below about IGeV/c^; 

171 4) a muon chamber system (MUG) made of 1000 m^ of 

172 Resistive Plate Chambers (RPG) arranged in 9 layers in 

173 the barrel and 8 layers in the end-caps and incorporated 

174 in the return iron yoke of the superconducting magnet. 

175 The position resolution is about 2 cm. 

176 The optimization of the event selection and the es- 

177 timate of physics backgrounds are performed through 

178 Monte Garlo (MG) simulations. The GEANT4-based 

179 simulation software BOOST [IJ] includes the geometric 

180 and material description of the BESIII detectors and the 

181 detector response and digitization models, as well as the 

182 tracking of the detector running conditions and perfor- 

183 mance. The production of the J/ip resonance is simu- 

184 lated by the MG event generator KKMG [H] , while the 

185 decays are generated by EVTGEN for known de- 

186 cay modes with branching ratios being set to PDG [l^l 

187 world average values, and by LUNDGHARM for the 

188 remaining unknown decays. The analysis is performed in 

189 the framework of the BESIII offline software system 



190 which takes care of the detector calibration, event recon- 

191 struction and data storage. 



192 II. EVENT SELECTION 

193 Signal J/ip — > Lupp events with lo — s- jtt'^ final states 

194 have the topology jjjpp. The event candidates are re- 

195 quired to have two well reconstructed charged tracks with 

196 net charge zero, and at least three photons. 

197 Gharged-particle tracks in the polar angle range 

198 I cosd\ < 0.93 are reconstructed from the MDG hits, only 

199 tracks in barrel region (| cos9\ < 0.8) are used to reduce 

200 systematic uncertainties in tracking and particle identi- 

201 fication. Tracks with their points of closest approach to 

202 the beamline within ±10 cm of the interaction point in 

203 the beam direction, and within 1 cm in the plane perpen- 

204 dicular to the beam are selected. TOE and dE/dx infor- 

205 mation are combined to determine particle identification 

206 confidence levels for tt, K and p{p) hypotheses; and the 

207 particle type with highest confidence level is assigned to 

208 each track. A proton and an anti-proton are required. 

209 To reduce the systematic error due to differences of the 

210 tracking efficiency at low momentum between data and 

211 MG, the momentum of the proton or anti-proton is fur- 

212 ther required to be larger than 300 MeV/c. 

213 Photon candidates are reconstructed by clustering sig- 

214 nals in EMG crystals. The photon candidates are re- 

215 quired to be in the barrel region (|cos0| < 0.8) of the 

216 EMG with at least 25MeV energy deposition, or in the 

217 end-caps region (0.86 < | cos 6*1 < 0.92) with at least 

218 50 MeV energy deposition, where 9 is the polar angle of 

219 the shower. Timing information from the EMG is used to 

220 suppress electronic noise and energy depositions that are 

221 unrelated to the event. To suppress showers generated by 

222 charged particles, the photon candidates are furthermore 

223 required to be separated by an angle larger than 10° and 

224 larger than 30° from the proton and anti-proton, respec- 

225 lively. 

226 A four-constraint (4G) energy-momentum conserving 

227 kinematic fit is performed to the jjjpp hypothesis. For 

228 events with more than three photon candidates, the com- 

229 bination with the minimum xic selected, and xlc 

230 is required. The tt*^ candidates are reconstructed from the 

231 two of the three selected photons with an invariant mass 

232 closest to the tt" mass, and |M^^ — A^^o| < 15MeV/c^ is 

233 required. 



234 III. BRANCHING FRACTION AND YIELD 

235 MEASUREMENTS 

236 

237 Figure 1 shows the 777° invariant mass spectrum for 

238 candidate J/V' — > jn'^pp events, where a distinctive uj 

239 signal is seen. An unbinned maximum likelihood fit is 

240 performed to the 771° invariant mass with the uj signal 
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M(y/) (GeV/c^) 

FIG. 1. 77r'' invariant mass distribution of J/ip —J- 'y-u'^pp 
candidates. The dashed hne is the signal shape which is 
parametrized by a Breit-Wigner function convoluted with 
the detector resolution described by the Novosibirsk func- 
tion; the dashed-dotted line is the background shape which 
is described by a second order Chebychev polynomial; and 
the solid line is the total contribution of the two compo- 
nents. The solid arrows indicate the uj signal region (0.753 < 
Mi^-K^) < 0.813 GeV/c^) and the two pairs of dashed arrows 
indicate the lo sidebands (0.663 < M{-fTr°) < 0.693 GeV/c^ 
and 0.873 < M(77r°) < 0.903 GeV/c^). 



241 parametrized by a Breit-Wigner function convoluted with 

242 the Novosibirsk function [30| which describes the detec- 

243 tor resolution. The background shape is described by a 

244 second-order Chebychev polynomial function. The mass 

245 and width of the w peak are fixed to the values published 

246 by the Particle Data Group (PDG) [l^], and the yield of 

247 the u! signal obtained from the fit is Nobs = 2670 ± 69. 

248 The branching fraction of J/ijj — s- cjpp is calculated 

249 according to : 



S( J/^ ^ upp) = ^^/^ ^ ^ 7.0) Vs(.o ^ X e 

. 

250 where iVobs is the number of signal events determined 

251 from the fit to the 77r'^ invariant mass; Nj/^, is the num- 

252 ber of J/ip events [l^l; B{uj -5- 777°) and B{tt° 77) 

253 are branching fractions of w — s- jn'^ and tt" — > 77, re- 

254 spectively, as from the PDG [l^l; and the detection efh- 

255 ciency Erec is (16.1 ± 1.7)% obtained from a MC sample 

256 for J/^ ^ ujpp events generated according to a phase- 

257 space distribution. The measured branching fraction is 

258 B{J/i} ujpp) = (9.0 ± 0.2 (stat.)) x IQ-". 



259 Candidate J upp events are selected with the 

260 mass window requirement 0.753 GeV/c^ < Af(77r") < 

261 0.813 GeV/c^; and the Dalitz plot of these events is shown 

262 in Fig. [21 There are no obvious structures in the Dalitz 

263 plot, though the distribution is different from the pure 

264 Lopp phase space distribution. The corresponding pp, 

265 Lop and up invariant-mass spectra are also presented in 

266 Fig. [2] The data points with error bars are from signal 

267 region and the hatched area are from the sideband region. 

268 the mass threshold is shown in Fig. [S] 

269 To obtain the number of J/^ — >■ uXijyp) ujpp 

270 events, an unbinned maximum likelihood fit is performed 

271 to the pp invariant mass around the mass threshold. In 

272 the fit, the spin-parity of X(pp) is assumed to be 0~, 

273 and the signal of X(pp) in the J/tp — > ujX{pp) — > ujpp 

274 decay is parametrized by an acceptance- weighted 5- wave 

275 Breit-Wigner function : 



BW{M) 



(Af2 - M2)2 4- AqV^ 



X ercc(M) . (2) 



276 Here, q is the momentum of the proton in the pp rest 

277 frame; k is the the momentum of the uj meson; i = 

278 is the relative orbital angular momentum; M is the in- 



variant mass of pp; Mq and F are the mass and width 
of the X{nn), respectively, which are taken from BESI- 
II results 3; Erec is the detection efficiency. The non- 
uj background is presented by a function of the form 
f{5) ^ N{S^/^ + aiS^/^ + a2S^/^) with 5 = Mpp - 2mp 
where rup is the proton mass. The normalization and 
shape parameters ai and 02 are determined by a simulta- 
neous fit to the M{pp) in uj signal region and uj sideband 
region 0.09GeV/c2 < |M(77r°) - 0.783| < 0.12GeV/c2. 
The non-resonant J/ijj ^ uipp events are also described 
by the function f{S), where the normalization and shape 
parameters are allowed to float. The fit results are shown 
in Fig. [3l and the number of X{pp) events is ± 1.6. 
A Bayesian approach p7| estimate the upper limit of 
B{J/ij ijX{pp) ujpp), and TVobs < 9 at 95% C. 
L. is determined by finding the value i^Dbs with 



= 0.95, 



(3) 



295 where Nobs is the number of signal events, and C is the 

296 value of the likelihood function with the Nobs value fixed 

297 in the fit. The upper limit on the product of branching 

298 fractions is calculated with 



B{J/iJj ujX{pp) ujpp) < 



Nj/^ X (1 



.) X B{uj -> 77r") X 6(7r0 —5- 77) x s^c 



(4) 
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FIG. 2. Dalitz plot and pp, uip, ujp invariant-mass spectra of J/ip — j- upp candidates. The data points with error bars are from 
signal region and the hatched areas are from the sideband region. 



2M where asys. is the total systematic uncertainty which will 

300 be described in the next section. The upper limit on the 

301 product of branching fractions is B{J/ip — )■ uiX{pp) — 5- 

302 upp) < 3.9 X 10"^ at the 95% C.L.. 

303 An alternative fit with a Breit-Wigner function includ- 

304 ing the Jiilich FSI 



f w „2L + lj,3 



305 for X{pp) is performed. Here, /fsi is the Jiilich FSI cor- 

306 rection factor The mass and width of Xipp) are 

307 taken from the previous BESIII PWA results The 

308 upper limit on the product of branching fractions is dc- 

309 tcrmincd to be BiJ/ip -> ioX{pp) ^ ujpp) < 3.7 x 10"^ 

310 at the 95% C.L.. 



311 IV. SYSTEMATIC UNCERTAINTIES 

312 Several sources of systematic uncertainties are con- 

313 sidcred in the measurement of the branching fractions. 

314 These include differences between data and the MC sim- 

315 ulation for the tracking algorithm, the PID, photon de- 

316 tection, the kinematic fit, as well as the fitting procedure, 

317 the branching fraction of the intermediate states and the 

318 total number of J/ijj events. 

319 The systematic uncertainties associated with the track- 

320 ing efficiency and PID efficiency have been studied with 

321 J/ijj ^ pp'!T^TT~ using a technique similar to that dis- 

322 cussed in Ref. [3l| . The difference of tracking efficiencies 

323 between data and MC simulation is 2% per charged track. 

324 The systematic uncertainty from PID is 2% per proton 

325 (anti-proton). 

326 The photon detection systematic uncertainty is studied 

327 by comparing the photon efficiency between MC simula- 

328 tion and the control sample J/ip — > pn. The relative 

329 efficiency difference is about 1% for each photon (32l. [33|. 
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FIG. 3. Near-threshold pp invariant-mass spectrum. The sig- 
nal J/ip LoX{pp) — >■ u)pp is described by an acceptance- 
weighted Breit-Wigner function, and and signal yield is con- 
sistent with zero. The dotted line is the shape of the signal 
which is normalized to five times the estimated upper limit. 
The dashed line is the non-resonant contribution described by 
the function f{5) and the dashed-dotted line is the non ujpp 
contribution which is estimated from uj sidebands. The sol- 
id line is the total contribution of the two components. The 
hatched area is from the sideband region. 

330 Here, 3% is taken as the systematic error for the efficien- 

331 cy of detecting three photons. The uncertainty due to tt" 

332 reconstruction efBciency is taken as 1% js^ [33j. 

333 To estimate the uncertainty associated with the kine- 

334 matic fit, selected samples oi J/ip -> E+S^ pTr^pn^ 

335 events are used. The kinematic fit efficiency is defined as 

336 the ratio between the signal yield of S+ with or without 

337 the kinematic fit. The difference of kinematic fit cffi- 

338 ciency between data and MC is 3%, and is taken as the 

339 systematic uncertainty caused by the kinematic fit. 

340 As described above, the yield of J/^p -> ujpp is de- 

341 rived from a fit to the invariant-mass spectrum of jtt'-' 

342 pairs. To evaluate the systematic uncertainty associ- 

343 ated with the fitting procedure, the following two as- 

344 pects are studied (i) Fitting region: In the nominal fit, 

345 the mass spectrum of jtt'^ is fitted in the range from 

346 0.663 GeV/c^ to 0.903 GcV/c^. Alternative fits within 

347 ranges 0.653 GeV/c^ to 0.913 GeV/c^ and 0.673 GeV/c^ 

348 to 0.893 GeV/c^ are performed, and the difference in 

349 the signal yield of 2% is taken as the systematic un- 

350 certainty associated with the fit interval, (ii) Back- 

351 ground shape: To estimate the uncertainty due to the 

352 background parametrization for the branching fraction 

353 B{J /tp — > i^pp), a first or third order instead of a second- 

354 order Chebychev polynomial is used in the fitting. The 

355 difference of 1.2% is used as an estimate of the systematic 

356 uncertainty. 

357 For the upper limit on the branching fraction B{J/'ip — ?• 



358 ll)X{pp) i^pp), the systematic uncertainty associat- 

359 ed with the fitting procedure is estimated by fixing 

360 the shape of the non-resonant contribution to a phase 

361 space MC simulation of J/?/' — ujpp, which is pre- 



90 F 




FIG. 4. Near-threshold pp invariant-mass spectrum. The sig- 
nal J/^ — ^ ujX(pp) ujpp is described by an acceptance- 
weighted Breit-Wigner function, and and signal yield is con- 
sistent with zero. The dashed line is the non-resonant contri- 
bution fixed to a phase space MC simulation of J /tj) — > uipp 
and the dashed-dotted line is the non upp contribution which 
is estimated from uj sidebands. The solid line is the total con- 
tribution of the two components. The hatched area is from a 
phase space MC simulation oi J /ip — >■ ujpp. 



362 sented by Figure. |31 enlarging/reducing the normaliza- 

363 tion of the non-cj contribution by 7% (the difference of 

364 the estimation of non-w background level between da- 

365 ta and inclusive MC); and varying the sideband region 

366 to 0.095 GeV/c2 < |M(77r°) - 0.783| < 0.115 GeV/c^ 

367 and 0.085 GeV/c2 < |M(77rO) - 0.783| < 0.125 GeV/c^. 

368 When fitting with or without the FSI effect, the signal 

369 yields for the alternative fits are lower or equal to the 

370 nominal fit, therefore the conservative upper limit from 

371 the fit without FSI correction is reported. 

372 Various distributions obtained with data and the 

373 phase-space MC sample have been compared and some 

374 discrepancies are observed. To determine the systemat- 

375 ic error on the detection efficiency associated with these 

376 discrepancies, an alternative detection efficiency is esti- 

377 mated by the re-weighting phase-space MC samples. The 

378 difference in detection efficiency compared to the nomi- 

379 nal one is 7% and taken as a systematic uncertainty. The 

380 number of J /ip events is determined from an inclusive 

381 analysis of J /ip hadronic events and an uncertainty of 

382 1.24% is associated to it [23|. The uncertainties due to 

383 the branching fractions of a; — )■ ^i:^ and tt*^ — 77 are 

384 taken from the PDG [13 ■ 



385 
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TABLE I. Summary of systematic 


uncertainties. 


'-' means the correspondinj 


5 systematic uncertainty is negligible. 






Upper limit of 


Upper limit of 


Source 


B(J/-rlj Ljpp) B(J/iP ^> ujX{pp) ^ 


ojpp) B{J/tl} '^X{pp) — >■ ujpp) with FSI 


Tracking 


4% 


4% 


4% 


PID 


4% 


4% 


4% 


Photon 


3% 


3% 


3% 


Kinematic Fit 


3% 


3% 


3% 


TT reconstruction 


1% 


1% 


1% 


Fitting region 


2% 




— 


Background Shape 


1% 






Branching fraction of intermediate state 


3% 


3% 


3% 


Total J/il> numbers 


1.24% 


1.24% 


1.24% 


MC Generator 


7% 






Total uncertainty 


10.3% 


7.8% 


7.8% 



386 V. SUMMARY 

387 In summary, using (225.3 ± 2.8) x 10^ J/'0 events col- 

388 lected with the BESIII detector, the decay of J/V' — > 

389 uipp in the decay mode uj ^i:'^ is studied. The 

390 branching fraction B{J/tjj — > t^pp) is measured to be 

391 (9.0 ± 0.2(stat.) ± 0.9(syst.)) x 10"^. No obvious en- 

392 hancement around the pp invariant-mass threshold is ob- 

393 served. At the 95% C.L., the upper limits on the prod- 

394 uct of branching fractions B{J /tp — > ioX{pp) — ujpp) are 

395 measured to be 3.7 x 10~^ and 3.9 x 10~® with and with- 

396 out accounting for the Jiilich FSI effect, respectively. As 

397 isospin for J/ip jpp and ujpp should both favor 1 = 

398 {1 = 1 should be suppressed in J/'ip jpp as in oth- 

399 er J/iJj radiative decays), the non-observation of X(pp) 
mo in Ljpp disfavors the pure FSI interpretation for the pp 

401 threshold enhancement in the decay J/ip ^ IPP- 
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